Programmable quantum state discriminators with simple programs 
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We describe a class of programmable devices that can discriminate between two quantum states. 
We consider two cases. In the first, both states are unknown. One copy of each of the unknown 
states is provided as input, or program, for the two program registers, and the data state, which is 
guaranteed to be prepared in one of the program states, is fed into the data register of the device. 
This device will then tell us, in an optimal way, which of the templates stored in the program 
registers the data state matches. In the second case, we know one of the states while the other is 
unknown. One copy of the unknown state is fed into the single program register, and the data state 
which is guaranteed to be prepared in either the program state or the known state, is fed into the 
data register. The device will then tell us, again optimally, whether the data state matches the 
template or is the known state. We determine two types of optimal devices. The first performs 
discrimination with minimum error, the second performs optimum unambiguous discrimination. In 
all cases we first treat the simpler problem of only one copy of the data state and then generalize 
the treatment to n copies. In comparison to other works we find that providing n > 1 copies of the 
data state yields higher success probabilities than providing n > 1 copies of the program states. 

PACS numbers: 03.67-a, 03.65.Ta, 42.50.-p 



I. INTRODUCTION 

Quantum state discrimination U is a basic tool for 
many tasks in quantum information and quantum com- 
munication. In the prototype problem a quantum pro- 
cessor generates a quantum system as its output which is 
in one of a set of known states but we do not know which 
and want to determine the actual state. If the possible 
states are not orthogonal this cannot be done with 100% 
probability of success since the cloning of quantum states 
is impossible. There are two basic strategies to accom- 
plish state discrimination. In the first, every time a mea- 
surement is performed we want to identify the state of 
the output with one of the possible states. Clearly, errors 
must be permitted and in the error minimizing strategy 
the optimum measurement is such that the probability 
of error is minimum. The case of discriminating with 
minimum error between two possible states was treated 
in the pioneering work by Helstrom 0. More recently, 
the interest was focused on the unambiguous discrimina- 
tion. In this strategy we are not permitted to make an 
erroneous identification of the state. The cost associated 
with this condition is that sometimes we fail to identify 
the state altogether. In the optimum strategy the prob- 
ability of failure is a minimum. The optimal value of the 
failure probability for two known and equally likely pure 
states was obtained by Ivanovic, Dieks and Peres (IDP 
bound, 0,0, 01). Later Jaeger and Shimony general- 
ized the IDP bound for arbitrary preparation probabili- 
ties of the states, i. e. for arbitrary prior probabilities of 
the two possible states. 

The actual state-distinguishing device for two known 



states depends on the two states, \ipi) and \ip2), k e. 
these two states are "hard wired" into the machine. An- 
other approach is to supply the information about the 
states to be distinguished as inputs, in particular as quan- 
tum inputs. That is, one encodes the information about 
the states one wants to distinguish into a quantum state, 
which is then a kind of quantum program, that is sent 
into the discriminator at the same time as the particle 
whose state is to be identified. The first such device was 
proposed by Dusek and Buzek 0. This device distin- 
guishes the two states cos(^/2)|0) ± sin(</>/2)|l), and the 
angle <f> is encoded into a one-qubit program state in a 
somewhat complicated way. The performance of this de- 
vice is good; it does not achieve the maximum possible 
success probability for all input states, but the average 
value of its success probability, averaged over the angle 
(j>, is greater than 90% of the optimal value. In a series 
of recent works Fiurasek et al. investigated a closely re- 
lated programmable device that can perform a von Neu- 
mann projective measurement in any basis, the basis be- 
ing specified by the program. Both deterministic and 
probabilistic approaches were explored @, and experi- 
mental versions of both the state discriminator and the 
projective measurement device were realized [sj- Sasaki 
et al. developed a related device, which they called a 
quantum matching machine 10] . Its input consists of 
K copies of two equatorial qubit states, which are called 
templates, and TV copies of another equatorial qubit state 
I/). The device determines to which of the two template 
states I/) is closest. This device does not employ the 
unambiguous discrimination strategy, but optimizes an 
average score that is related to the fidelity of the tern- 
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plate states and |/). Programmable quantum devices to 
accomplish other tasks have been explored by a number 

of authors [13-123 ■ 

Recently two of us proposed an approach to a pro- 
grammable state discriminating machine in which the 
program is related in a simple way to the states \tj)-i) 
and \1jj2) that one is trying to distinguish |18| . A mo- 
tivation for this problem is that the program state may 
be the result of a previous set of operations in a quan- 
tum information processing device, and it would be easier 
to produce a state in which the information about l^i) 
and 1^2) is encoded in a simple way than one in which 
it is encoded in a more complicated way. The program 
is the most elementary possible, it consists of copies of 
the states one is trying to distinguish. The device then 
performs optimally with the given program states or, in 
other words, it optimally identifies the data state with 
one of the two unknown program states, or reference 
states, respectively. Despite the complete lack of classical 
information about the reference states, the identification 
is still possible due to symmetry properties that are in- 
trinsically quantum mechanical, and are similar to those 
first employed by Barnett et al. J2~3| for the purpose of 
comparing unknown states. 

The original results in |l8l ] were for unambiguous dis- 
crimination and for qubit data and program states. They 
have recently been extended to qudits by Hayashi et al., 
both for optimum unambiguous discrimination |19| and 
for minimum error discrimination |2(ij |. Their investiga- 
tions are restricted to equal prior probabilities, where the 
data state equally likely matches each one of the program 
states, but they also dealt with the case in which an ar- 
bitrary number of copies is provided for each of the two 
program states. 

In the present paper we generalize the programmable 
state discriminator introduced in [Tsl ] in several other 
directions and develop a comparative study of pro- 
grammable state discriminators based on the two mea- 
surement strategies of minimum-error discrimination and 
optimum unambiguous discrimination. For this pur- 
pose in Sec. II we first reformulate the problem of pro- 
grammable state discriminators as a problem of discrim- 
ination between two mixed quantum states. Sec. Ill 
is devoted to the case that both of the pure states to 
be discriminated are unknown so we need a reference 
state for each. In Part A we treat the error minimizing 
version of the programmable state discriminator, con- 
sidering both a joint measurement on all three qubits, 
and also a measurement prescription that is restricted 
to two-qubit measurements only. In Part B we rederive 
the results of Ref. 18] for the unambiguous version of 
the programmable state discriminator partly for compar- 
ison's sake but also using the consistent approach based 
on the equivalent mixed state discrimination problem. It 
should be noted, in this context, that the results of Rcf. 
|18| were obtained in a somewhat ad hoc manner and 
the current approach gives a solid foundation to those 
results. We also compare the optimal probabilities ob- 



tained in Parts A and B for the programmable state dis- 
criminators based on the two possible strategies. In Sec. 
IV we fill another gap and show how to construct devices 
that can optimally discriminate between one known and 
one unknown state using both minimum-error and op- 
timum unambiguous strategies. That is, we know what 
\tpi) is, but do not know \ip2)- Then we need a refer- 
ence state only for the unknown state, which constitutes 
the program in this case. We can say that this line of 
investigation characterizes the quality of the source that 
produces the states to be discriminated, or the quality 
of our knowledge about the source, respectively. If both 
possible states are known (the original IDP and Helstrom 
problem) there is no need for a program, the states are 
hard wired into the optimal device. If one of the states 
is known we need a program for the unknown state while 
the other is hard wired into the device and if both states 
are unknown we need a program for both. 

We also take look at another aspect of the problem. 
Namely, besides investigating the effect of the source 
quality on the optimal performance of this family of state 
discriminating devices, we also investigate the effect of 
the resources on the performance of these devices. Sup- 
pose that instead of one copy of the state to be discrimi- 
nated we are given n copies, but we still only possess one 
copy each of the unknown state(s) as the program state 
(or none for two known states). In Sec. V we therefore 
generalize the two unknown qubit scenarios of Sec. Ill for 
the case when n copies of the input state, and one copy of 
the program states, are provided. In Sec. VI we provide 
a similar generalization of the one unknown qubit cases 
treated in Sec. IV. In each of these cases we determine 
the optimal measurement strategy both for minimum er- 
ror and unambiguous discrimination of the data state. 
The devices that accomplish this are programmable, in 
the first case the program consists of two qubits, one in 
\tpx) and one in IV'2), while in the second case the pro- 
gram consists of a single qubit in the state \ip2)- Note 
that in all cases, the program is extremely simple. It is 
what could be called a "quantum list", a set of qubits, 
one in each of the states to be discriminated, or one each 
in some subset of the states to be discriminated. In Sec. 
VII we conclude with a brief discussion of how these re- 
sults can be used to characterize the preparation quality 
(source quality) and to quantify the available resources. 



II. DISCRIMINATION OF UNKNOWN STATES 
AND ITS CONNECTION TO THE 
DISCRIMINATION OF MIXED STATES 

Let us begin by briefly reviewing the problem that was 
originally addressed in 18] . We consider a system of 
three qubits, labeled A, B, and C, and assume that the 
qubit A is prepared in the state and the qubit C is 
prepared in the state {1^2) ■ Qubit B is guaranteed to be 
prepared in either or | -02 ) , with a probability 771 of 
being in \ipi) and a probability 772 = 1 — ?/i of being in 
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| -02 ) - The states \tp%) or \ip 2 ) are different and unknown. 
Our task is to find whether the state of qubit B is \tpi) 
or | "02} - One way of looking at this problem is to view 
the qubits A and C as a program. They are sent into the 
program register of some device, called a programmable 
state discriminator, and the third, unknown qubit is sent 
into the data register of this device. The device then tells 
us, with an optimal probability of success, which one of 
the two program states the unknown state of the qubit 
in the data register corresponds to. We can consider 
this problem as a task in measurement optimization. We 
want to find an optimal measurement strategy that, with 
a maximum probability of success, tells us which one of 
the two program states, stored in the program register, 
matches the unknown state, stored in the data register. 
In [Tsj only unambiguous discrimination was treated, in 
which the measurement is allowed to return an inconclu- 
sive result but never an erroneous one. Here we want to 
investigate the measurement strategy of minimum-error 
discrimination, as well. In general, we want to determine 
the best possible measurement for identifying the state 
of the qubit B. Our task is then reduced to the following 
measurement optimization problem. One has two input 
states 

|*l) = \^i)a\^x)b\^2)c , 

1*2) = \ih)A\ih) B \ih)c , (2.1) 

where the subscripts A and C refer to the program reg- 
isters (A contains and C contains ^2) )> and the 
subscript B refers to the data register. Our goal is to 
optimally distinguish between these inputs, with respect 
to some reasonable criteria, keeping in mind that one 
has no knowledge of \ipi) and \tp 2 ) beyond their a priori 
probabilities. 

Assuming the states j^i) and \4>2) to be completely 
unknown, we have to find the measurement strategy that 
is optimal on average. Thus, we have to take the average 
of the input with respect to all possible qubit states. The 
problem is then equivalent to distinguishing between two 
mixed states, given by the density operators 

Pi={|*i)(*i|} av , (2.2) 

P2 = {|*2>(*2|} av , (2.3) 

that occur with the prior probabilities 771 and 772, respec- 
tively. Any state of a particular qubit (A, B or C) can 
be represented using the Bloch parametrization given by 
\ifii) = 008(^/2)10) +e 1 *sin(0 i /2)|l) (i=l,2), with |0) 
and |1) denoting an arbitrary set of orthonormal basis 
states. Here 9 and <f> are the polar and azimuthal an- 
gle on the Bloch sphere. After performing the averaging 
with respect to all possible values of and <p we arrive at 

Pi = \p S a V b®Ic, (2-4) 
92 = \l A ®PT ■ (2-5) 



where P^b = J2i=i \ u i)ac Ac{ui\ and P B V ™ = 
Si=i \ u i)BC BC'( u i\ are the projectors onto the symmet- 
ric subspaces of the corresponding qubits, AB and BC, 
respectively. Here we used the two-qubit basis states 

I \ m\ ln\ l \ \0)a\1)b + \1)a\0)b ,„ 
\ui)ab = |0)^|0)b, \u 2 )ab = -j= , (2.6) 

I \ h\ h\ |-\ |0)a|1)b-|1)a|0)b , 97 x 
\u 3 )ab = \1)a\1)b, \u 2 )ab = ^= ,(2.7) 

and the analogous expressions for the qubit combination 
BC. Due to the symmetry of the state p\ with respect 
to interchanging qubits A and B the antisymmetric state 
\uijAB does not enter the expression for the density oper- 
ator. Eqs. (|2.4|l and 12.5fl reduce our state identification 
problem to the problem of discriminating between these 
two mixed states. 

When one of the two states that we want to distinguish 
is known, we arrive at a simpler variant of the discrimi- 
nation problem. We do not need to provide a template 
for the known state, and one of the program registers, 
say A, can be eliminated from the problem. It is con- 
venient to define the single-qubit basis states in such a 
way that the known pure state serves as one of the basis 
states, denoted by |0), so = |0). We then have to 
distinguish two cases, the qubit B is either in the state 
|0)b, occurring with the prior probability r]i, or it is in 
the unknown state of the qubit C, occurring with the 
prior probability 772. These two cases correspond to the 
density operators 

Pi = \o)B B {o\®MccM} m = l\o)BB{o\®ic 

= 2 {\ u ^)bc bc(ui\ + \v 2 )bc bc(v2\} , (2.8) 

P'2 = Mb \^)c bW c(^|} av = \P B V c- (2-9) 
Here 

\v 2 )ab = \0}a\1)b- (2.10) 

In addition, we introduce 

\v 2 )ab = \1)a\0)b, (2.11) 

together with the analogous expression for the qubit com- 
bination BC. {|u2)> 1^2)} and{|u2), ^2)} form alternative 
bases for the subspace with exactly one qubit in the state 
|1). They will prove useful later when we consider the 
various discrimination scenarios in the following sections. 

After these preliminary considerations we are now in 
a position to investigate different possible measurements 
for identifying the state of the qubit B, i. e. for distin- 
guishing between the two density operators given by Eqs. 
(|2~3|) and i(2~5|) or, alternatively, by Eqs. l|2~%j) and (|2~3|) . 
Before doing so, we briefly recall the underlying theoret- 
ical concepts for treating the strategies of discriminating 
two mixed states. Any measurement suitable for distin- 
guishing between the mixed states p± and p 2 , occurring 
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with the prior probabilities r/i and 772 = 1 — 771, respec- 
tively, can be formally described with the help of three 
positive detection operators Ilo, 111 and II2, whose sum 
is the identity, 



n + Ei + n 2 = 1. 



(2.12) 



These operators are defined in such a way that for j = 1,2 
Tr(pllj) is the probability to infer from the measurement 
that the system is in the state pj if it has been pre- 
pared in a state p, while Tr(pE ) is the probability that 
the measurement result is inconclusive, i. e. that the 
measurement fails to give a definite answer. When all 
detection operators are projectors, the measurement is 
a von Neumann measurement, otherwise it is a gener- 
alized measurement based on a positive operator-valued 
measure (POVM). Once the detection operators of a gen- 
eralized measurement have been found, Neumark's the- 
orem guarantees that schemes for actually realizing the 
measurement can be devised by determining suitable pro- 
jections in an enlarged Hilbert space that results from 
appending an ancilla to the original system [24I |25| . 

The above POVM is appropriate for unambiguous 
state discrimination. For minimum-error discrimination 
inconclusive results do not occur, so that 



n = 0, 



(2.13) 



and we require that the probability of errors in the dis- 
crimination procedure is a minimum. For two mixed 
states this problem was originally solved by Helstrom [2j . 
The error probability is always larger than zero unless the 
states to be distinguished are orthogonal, and it can be 
expressed as 

Perr = r?lTr(piE 2 ) + 772Tr(/9 2 Ei) 

= m + Tr[(772p 2 - niPlW, (2-14) 

where in the second line Eqs. (|2.12Jl and l|2.13|l have been 
used, as well as the relation 772 = 1 — 771. After introducing 
the operator 



(2.15) 



it is obvious that the minimum of the error probability is 
obtained when n! is the projector onto those eigenstates 
\4>k) of A that belong to negative eigenvalues The 
optimum detection operators therefore read 



n opt 



k<k 



k > kn 



where A& < for 1 < k < fco and A^ > for k > 
ko- Clearly, these two operators are projections, and the 
optimal minimum-error measurement for discriminating 
between two quantum states is, therefore, always a von 
Neumann measurement. The resulting minimum error 
probability P™™ = Pe is given in Q by 



p E = 2 (1 - Tr|7? 2 p2 - mPi\) = 2 



!(.-!>,) 



In optimum unambiguous discrimination which is the 
other frequently used strategy errors are not allowed to 
occur. This requirement is equivalent to 

P11I2 = p 2 ni = 0, (2.18) 

(see, for example lj). In the optimum measurement 
scheme the failure probability, i. e. the probability for 
getting an inconclusive outcome, is minimized, taking 
into account the constraint that the eigenvalues of the 
operator no = I — Hi — H2 are non-negative. The fail- 
ure probability is always nonzero unless the states to be 
discriminated are orthogonal. It can be expressed as 

Qfaii = 77iTr(piE ) + 772Tr(/j 2 Eo) 

= 1 - 77iTr(/?iEi) - 77 2 Ti'( / 02E 2 ) 

1 Psucc: (2-19) 

where Eqs. (|2.12(l and l|2.18(l have been used, and where 
we also introduced the success probability P SUC c of the 
measurement. Optimum unambiguous discrimination 
between two mixed states is an issue of ongoing theo- 
retical research [SB, E3, H^, [H^ ES| In contrast 
to minimum-error discrimination, there does not exist a 
compact formula expressing the minimum probability of 
inconclusive results, i. e. the minimum failure probabil- 
ity, for unambiguously discriminating two mixed states 
that are completely arbitrary. However, analytical solu- 
tions can be obtained for certain special classes of density 
operators, including the cases that are of interest for this 
paper. 



III. TWO-QUBIT PROGRAM, SINGLE COPY 
OF THE DATA STATE 

A. Minimum-error discrimination strategy 

1. Joint measurement on all three qubits 

We begin by investigating the measurement that dis- 
criminates, with minimum probability of error, between 
the density operators given by Eqs. Ij2.4|l and Q2.5JI . For 
this purpose we define the orthonormal basis states in 
the eight-dimensional Hilbert space spanned by the three 
qubits as 



|1) = |0)a|0) b |0) c , 

|3) = |0)a|1)b|0)c, 

|5) = \0) A \1) B \1) C , 

|7) = \1)a\1)b\0)c, 



|2> = |0)a|0) b |1>c, 
|4) = |l) A |0) B |0)c, 
|6) = |1) A |0) S |1) C , 
|8) = |1) A |1) S |1) C . (3.1) 



The numbering of the states is essentially the binary 
number formed by the bit values on the right-hand side 
shifted by one. Note, however, that in the case of |4) and 
1 5) the order is reversed. Expanding the expressions for 
pi and P2 in this basis and introducing the notations, 



(2.17) 



In) 



13) + |4) 
V2 • 



15) + |6) 
V2 



(3.2) 
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|2> + |3> 



s 2 = 



|7> 



V2 ' ' ' 
we obtain the spectral representations 

2 

1 ,i 

Pi 



(3.3) 



P2 



1=1 

= \si)(si\ + |1)<1| + |4)<4| + |5)<5| + |8)(8|).(3.5) 



When we express the operator A = T] 2 p2 — ViPi with the 
help of the basis states given by Eq. (|3.1|) , we arrive at an 
eight-dimensional square matrix which is block-diagonal 
if the columns and rows are numbered according to the 
numbering of the basis states. It can be written as 



A 



/ Li 

L 3 

L 3 

V Li 



(3.6) 



where L\ — {r\i — T)i)/6 and 



L3 =12 



'/2 



m 
o 



o 

-m 



m - m 
-m 2r ?2 - vi 



Since the eigenvalues of L3 are (772 — r/i) /6 and 



m - m 



± V (% - f?i) 2 + 3r?i^2 



(3.7) 



(3- 



the complete set of eigenvalues of the operator A is given 
by 



Ai — A2 — A_, A3 — A4 — A+, 

m - m 



A fc = 



G 



(5 < k < 



(3.9) 



The corresponding eigenstates, \4>k), are found to be 
\<Pr) = |1), |0 8 ) = |8), 

1 , . _ o~|2) + |3) + 6~|4) 



103) 
105) 



v/l + (a-) 2 + (&-) 2 
q+|2) + |3)+b+|4) 
y/1 + (a+) 2 + (6+) 2 
|2) + |3) + |4) 
v/3 



(3.10) 



where 



a ± = ±\/{r)i - r] 2 ) 2 + 3771772 - rji 
b ± = T\/ (m ~ V2) 2 + 3771772 - ?72 



(3.11) 



The eigenstates |02), 104), 106) follow from replacing the 
ordered set {|2>, |3), |4)} by the ordere d set {|5>, |6>, |7)} 
in the expressions given by Eq. (|3.1U|) . By inserting the 
eigenvalues of A into Eq. (|2.17|l we find after a little 



algebra that the minimum probability of error can be 
written in the following compact way, 



Pe = Vn 



1 



1 



(3.12) 

where r] max (rj m i n ) is the larger (smaller) of 771 and 779. 
This result is in agreement with the one derived in |2Qj 
for the special case that r] ma x — ilmin = 1/2. 

The above expression lends itself to a transparent in- 
terpretation. The error probability P e would be 77 m i„ if 
we did not perform any measurement at all but would 
simply guess, always choosing the state whose a priori 
probability is larger. The factor in the bracket, multi- 
plying r\ m ini is the improvement due to the optimized 
measurement. It is a slowly varying function of the prior 
probabilities, its value lying between 0.71 and 0.75 in the 
entire < 77 mi „ < 1/2 interval. To be specific, let us 
assume that the qubits are labeled in such a way that 771 
is the smaller of the two prior probabilities, i. e. that 
771 < 0.5. In this case Ai and A2 are the only negative 
eigenvalues, and the optimum detection operators II° pt 
and Il2 Pt for minimum-error identification take the form 
given by Eq. (|2.1(i|l with fc = 3. This means that the 
qubit B is inferred to be in the state of qubit A when a 
projection onto the subspace spanned by the eigenstates 
of II° pt is successful, and after successful projection onto 
the complementary subspace it is inferred to be in the 
state of qubit C. 

From the structure of the eigenstates \4>k) determining 
the optimum detection operators, and from the definition 
of the basis states, given by Eq. I|3.1|l . it is obvious that 
the smallest possible error probability, Pg, can only be 
obtained by performing a joint measurement on all three 
qubits simultaneously. The question therefore naturally 
arises as to what is the smallest value of the error prob- 
ability achievable under the restriction that only joint 
measurements on two qubits are allowed. In the following 
we study this problem. This situation is worth examining 
for two reasons. First, two-qubit measurements are easier 
to perform than three-qubit ones. Second, by comparing 
the results of the two- and three-qubit measurements, we 
see how the additional quantum information contained 
in the third qubit affects the result. 



2. Restriction to two-qubit-measurements 

First let us assume that the qubit C is not accessible, 
but that we are able to perform a joint measurement 
on the qubits A and B. This would be the case, for 
example, if a copy of only one of the two states we are 
trying to distinguish is provided. Starting again from 
Eqs. (|2.4|) and (|2.5|) . the problem of identifying the state 
of B is then equivalent to discriminating between the two 
reduced density operators, 



Tic Pi 



1 psym 

3 AB 



(3.13) 



6 



Tr c p 2 



1 



Ib 



(3.14) 



with the prior probabilities 771 and 772, respectively. 

The operator A = rj 2 Trcp 2 — ^lTrcpi, relevant for 
minimum-error discrimination, is now given by 



u 2 \ 



A = - y) XI \ u i)ABAB(Ui\ + y^ABAB- 

i=l 

(3.15) 

When ?7i < 3772/4 (or 771 < 3/7) all four eigenvalues of A 
are positive, and from Eq. (|2.16(l the optimum detection 



operators are obtained as TEf* — I and II° pt = 0. Hence 



the minimum error probability is achieved by guessing 
that the quantum system is always in the state that is 
more probable, in this case \i^ 2 ), without performing any 
measurement at all. This is a special situation, described 
earlier |33j. which has been observed in connection with 
a different problem of two-qubit-discrimination [34l | . 
On the other hand, for 771 > 3i] 2 /4: (or 771 > 3/7) 

\u 2 )abab{u 2 \ and II° pt = 



we readily find that II^ 15 * = 



J2i=i \ u i) ab ab{ui\, i. e. that the error probability is 
smallest when the qubit B is guessed to be in the state 
of qubit A after a successful projection onto the symmet- 
ric subspace of qubits A and B, while there is no guessing 
involved after a successful projection onto the antisym- 
metric subspace of qubits A and B. It is then known 
with certainty to be not in the state of qubit A. 

The results for the minimum error probability, follow- 
ing from Eqs. (|3.15|l and l|2.17[l . can be summarized as 



pAB _ 



l if Vi < f 

(1 — 771) otherwise. 



(3.16) 



Similarly, a joint measurement on the qubits B and C 
yields the minimum error probability 



pBC 



1 — rji otherwise. 



(3.17) 



Fig. 1 also reveals that by performing the optimal two- 
qubit measurement an error probability can be achieved 
that is almost as low as the absolute minimum error prob- 
ability, Pe, given by Eq. I|3.12[l where the latter can only 
be reached by a joint measurement on all three qubits. 
Even when the advantage of the three-qubit measure- 
ment is largest, which happens for equal prior probabili- 
ties, 771 = 772 = 1/2, the difference in the respective min- 
imum error probabilities for state identification is only 
marginal, 



P, 



AB 



P. 



BC 



0.375, Pi 



E — 7T 



4\/3 



= 0.356. (3.18) 




0.2 0.4 0.6 0.8 1 

FIG. 1: Comparison of the optimum performances of the er- 
ror minimizing strategy and the unambiguous discrimination 
strategy for the discrimination of two unknown pure states. 
The minimum error probabilities Pe resulting from a three- 
qubit measurement (full line) and P^ B , Pe C resulting from 
two-qubit measurements (dash-dotted and dash-double dot- 
ted line, respectively) are compared to the failure probabilities 



for unambiguous discrimination Q F B , Qp V (dotted lines) and 
qPOVm (hashed ii ne ). The error and failure probabilities are 
plotted vs. the prior probability 7/1 that the state of the data 
qubit B matches the state of the program qubit A. 



B. Optimum unambiguous discrimination strategy 

The optimum measurement for unambiguously iden- 
tifying the state of the data qubit B was found in [l^ 
using a method that relied on a special Ansatz for the 
detection operators, justified by the symmetry properties 
of the inputs. For completeness, here we reconsider the 
problem in the framework of the optimum unambiguous 
discrimination of two mixed states. 

In the following we apply the method developed in [3(j 
which, in turn, is a special case of the more general ap- 
proach in [32^. Starting from Eqs. (|3.4|l and 13.5|) . we 
denote the projectors onto the supports of p\ and p 2 by 
Pi and P 2 , respectively. The eigenstates of the operators 
I — P 2 and I — Pi are easily found to be 



01 



|2) - |3) 

|3) - |4) 
V2 • 



\°2) = 

N = 



17) 



V2 
15} - |6) 



V2 



(3.19) 
(3.20) 



Clearly, p 2 \a,i) = and p\\bi) = for i = 1,2. The 
most general Ansatz for the detection operators, satis- 
fying Iiip 2 = II2P1 = as required for unambiguous 
discrimination, therefore reads 



In the next paragraph we compare the minimum prob- 
abilities of error with the minimum probability of fail- 
ure arising in the other important measurement strategy, 
that of unambiguous discrimination. 



IF = J2 <xii\"i)("il n 2 = E 0io\h}(bj\. (3.21) 
From Eq. I|2.19|l we readily find that these two detection 
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operators give rise to the failure probability 



Qf 



ail 



= 1 -oBw 



(3.22) 



Note that due to the structure of the two given density 
operators the failure probability does not depend on the 
off-diagonal elements of the detection operators given by 
Eq. i|3.21[) . a property that is common to all problems 
of optimum unambiguous discrimination of two mixed 
states that have been explicitly solved so far |3(], 
l32l |. We are therefore free to choose = /% = for 
i ^ j, a choice that guarantees that IIo is positive for the 
largest possible values of an and pa (i = 1, 2), i. e. that 
Q fan can be made as small as possible. 

When we represent the operator IIo i n the basis defined 
in Eqs. (|3.1[1 . we again arrive at a block-diagonal eight 
by eight matrix, similar to Eq. (|3.6|l . given by 



n n = 




o o 

o o 

M(p 22 ,a 22 ) 

1 



(3.23) 



Here we introduced the abbreviation 

2 - a 



M{a,b) = \ 




(3.24) 



The eigenvalues of M(a, b) are found to be 1 and 
1 



H±{a,b) = 1 - - [a + b± y/{a- b) 2 + abj , (3.25) 

where obviously p±(a,b) — (i±(b,a). All eigenvalues of 
IIo are nonnegative provided that (j,+ (aa, Pa) is non- 
negative for i = 1,2 which holds true when j3a < 
(4 — 4au)/(4 — 3au). Hence in order to make Qf a u as 
small as possible, while keeping Ho a positive operator, 
we chose the equality sign and put 



Pi: 



4 - 4aj. 
4 - 3a f 



1,2). 



(3.26) 



After inserting these expressions into Eq. I|3.22|l the re- 
sulting function of an and a 22 has to be minimized, tak- 
ing into account that < an < 1, which follows from the 
fact that Tr(pHi) describes a probability for any density 
operator p. In accordance with the optimization problem 
solved in 0], we nn d that the failure probability takes 
its smallest possible value when an = a 22 = a, where 







I'2 



if vi < \m 
£) if km < m < 4^2 

if Vi > 4?7 2 . 



(3.27) 



Using Eqs. Ij3.21|) and (|3.1|) we arrive at the optimum 
detection operators 



n opt 



yrOpt 

11 2 



a I A ® \u 2 )bc (u 2 \bc, (3.28) 
4 — 4a 

-. 5- \u 2 )ab {u 2 \ab ® Ic, (3-29) 

4 — 3a 



where the value of a in the different parameter regions 
for r/i and t] 2 is given by Eq. (|3.27|) . Clearly, in the first 
parameter region n° pt = 0, while describes a pro- 
jection onto the antisymmetric two-qubit state \u 2 )ab- 
Similarly, in the third parameter region a projection onto 
the antisymmetric state \u 2 )bc has to be performed for 
optimum unambiguous discrimination. The failure prob- 
abilities resulting from these two von Neumann measure- 
ments are [13 



Q A F B 



m 
4 



IF 



1 



\ (3-30) 



In the intermediate region of the prior probabilities the 
optimum measurement is a generalized measurement, 
yielding the failure probability |isj 



n POVM 



3 (5<m<f), (3-31) 



where we took into account that r]2 = 1 — Tji . 

As seen in Fig. 1, the minimum failure probability is 
always at least twice as large as the minimum error prob- 
ability Pe for identifying the qubit state. This agrees 
with the general relation between the failure probabil- 
ity of optimal unambiguous discrimination and the er- 
ror probability of minimum-error discrimination of two 
mixed states that was derived in (35|. Fig. 1 also 
shows that the advantage of performing a generalized 
measurement, as compared to the best projective two- 
qubit measurement, is small. For 7/1 = r\% = 0.5, where 
QjP = QBC = 7 / 8 = 0j875 the failure probability is 

only reduced to the value Q p F OVM = 5/6 = 0.833. Of 
course, the surprise is not the high value of the failure 
probability but that the success probability is finite for 
the discrimination of completely unknown states. As we 
shall see in the next section, these relatively large fail- 
ure probabilities are somewhat reduced when one of the 
reference states is known. 



IV. ONE-QUBIT PROGRAM, SINGLE COPY 
OF THE DATA STATE 

A. Minimum-error discrimination between one 
known and one unknown state 



Now we treat the simplified case in which we want to 
decide whether the qubit B is in the known state |0)s, or 
whether it is in the unknown state of the program qubit 
stored in register C. We then have to distinguish between 
the density operators p\ and p' 2 given by Eqs. I|2.8|l and 
(|2.9|) that refer to the first and second alternative, respec- 
tively, and that occur with the prior probabilities 771 and 
772 • The subsequent treatment proceeds along exactly the 
same lines that we followed in the previous sections. 

For minimum-error identification we have to determine 
the eigenvalues and eigenstates of the 
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operator A' = r/2/o' 2 — f]\p'\- It is easy to obtain the 
spectral representation 

A ' = (f _ f) |ui>(Ml1 + f |U3>(W3 ' + ? ^i^x^i. 

(4-1) 

where A± = yj ( 2ry 2 — 3t7i ± \/ 4rj| + 97/2 J and 



According to Eq. (|2.1fij) it is reached with the help of the 
detection operators 

!\f-)Bc( i p-\Bc if m < I 

K)bc(wi|bc + |v-)-bc(V-|bc otherwise, 

(4.4) 

and n° pt = I BC - II° pt , where we have made use of 

the identity I B c = \ui){u\\ + 1^3) (w 3 | + J2i=± \<Pi){<Pi\- 
Clearly, the measurement that identifies the state of the 
qubit B with the smallest possible error is a joint projec- 
tion measurement on the qubits B and C . 

To close this section we briefly investigate the case 
that we can only perform a measurement on the qubit 
B alone, which means that our identification problem 
amounts to discriminating between the density opera- 
tor p'l = |0)bb(0| and the uniformly mixed state p' 2 ' — 
TrcP 2 = \^b- For minimum-error identification we 
consider the eigenvalues and eigenstates of the operator 
A" = (772/2 - m ) |0) SB (0| + 773/2 |1)bb(1|, obtaining in 
the standard way the minimum error probability 

P" = { f if Vl ~ 5 (4 5) 

E \ |(1 — ?7i) otherwise. 

If 771 < 1/3 the error probability is smallest when the 
qubit B is always guessed to be in the unknown state 
of the qubit C, while otherwise a projection measure- 
ment characterized by the operators II° pt = |0)_bb(0| 
and II 2 pt = |1)bb(1| has to be performed in order to 
minimize the error. 

As can be seen from Fig. 2, for 771 < 2/5 the optimum 
two-qubit measurement leads to a substantial reduction 
of the minimum error probability in comparison to the 
optimum single-qubit measurement. The difference is 
largest for 7/1 = 1/3, where P' E = 0.22 but P£ = 0.33. 



By making use of Eq. I|2.I 7(1 we find that the minimum 
error probability for identifying the state of the qubit B 
is given by 



(4.3) 

I 



1 




FIG. 2: Failure and error probabilities for the various strate- 
gies of discriminating between one known and one unknown 
state vs. the prior probability 771 . The minimum error proba- 
bilities P' E resulting from a joint measurement on the qubits 
B and C (full line) and from a single-qubit measurement, P'e 
(dash-dotted line) are compared to the failure probabilities 
for unambiguous identification Q'p , Q'p C (dotted lines) and 
Q'i OVM (dashed line). 



B. Optimum unambiguous discrimination between 
one known and one unknown state 

Finally we want to compare the minimum error prob- 
abilities with the minimum failure probability in unam- 
biguously identifying the qubit state. For this purpose 
we again use the method described in [3(j and [32 • By 
taking a reduction theorem |26| into account, from Eqs. 
(|2.8|) and (|2.9|l it follows that the most general Ansatz 
for the detection operators can be written as 

IIi = a|tia>(«a|, H 2 = /3\v 2 ){v 2 \ + |« 3 )(« 3 |, (4.6) 



\u 3 )~c ± \u 2 ) 2772 ± V^'j + 977? 
\f±) = ; c± = . (4.2) 



P'e = Tin 



max T}min) ^Vmin'H'mo, ? 



9 



where |«2) = (\ u 2) — |w2))/v2- Here again the subscript 
BC has been dropped. Clearly, IIip 2 = TI2P1 = as re- 
quired for unambiguous discrimination. As follows from 
Eq. H2.19fl . these detection operators yield the failure 
probability 



V2 , 



W = l-^«-^(2 + /3) 



(4.7) 



which has to be minimized under the constraint that 
IIo = 7 — IIx — II2 is a positive operator. For Ilo we 
obtain the expression 



Ho = I - a\u 2 ){u 2 \ - \u 3 )(u 3 \ 

\u 2 ))((u 2 \ ~ (u 2 \ 



0(\ \ 



(4.8) 



The eigenvalues of IIo are [ii = 1, jjb-z = and fi± — 
(2 — (3 — a ± \J cP- + (3 2 )/2, and they all are non- negative 
provided that (3 < (2 — 2a)/(2 — a). In order to minimize 
Qfaii while keeping Hq a positive operator we therefore 
choose 



2a 



(4.9) 



Upon substituting (3 into Eq. 1|4.7[) and determining the 
smallest value of the resulting function of a, taking into 
account that < a < 1, we find that the minimum failure 
probability is obtained when 



a 



if 37/1 < Tj 2 

if r/2 < 3?7i < 4?7 2 
if 3771 > 4?72. 



(4.10) 



Using Eqs. (|4.1UI) and (|4.9|l in Eq. (|4.6[) yields the ex- 
plicit expressions for the optimum detection operators. 

If 3?7i < 772, which implies that r\\ < 1/4, we have 
II° pt = and n° pt = |1) BB (1| O I c , i- e. the opti- 
mum measurement is a projection measurement on the 
qubit B alone. On the other hand, for Zrji > Ari2, i.e. 
Vi ^ 4/7, the optimum measurement is a joint pro- 
jection measurement on the qubits B and C, where 
n° pt = \u2)bcbc(u2\ and H° 2 pt = ® |l)cc(l|- 

The failure probability of these two von Neumann mea- 
surements is given by 



IB 

F 



Q 



= 1 
= 1 



1 - 771 

2 ' 
m _ m 
4 3 



m 
12' 



(4.11) 



In the intermediate parameter region the optimum mea- 
surement is a generalized measurement, yielding the fail- 
ure probability 



Q 



iPOVM 
F 



6 



The benefit of performing the generalized measure- 
ment is only marginal, as evidenced by Fig. 2. In fact, 



the reduction of the failure probability compared to the 
best of the two types of von Neumann-measurements is 
largest for 771 = 0.4, where Q'P = Q'f c = 0.7, while 
the generalized measurement yields the failure probabil- 
ity Q'p OVM — 0.683. In agreement with the general re- 
lation derived in |35| , the latter value is more than twice 
as large as the minimum error probability P' E which for 
771 = 0.4 takes its maximum value 0.258. 

A comparison of Figs. 1 and 2 reveals that, as ex- 
pected, the minimum probabilities of error and failure 
are smaller when one reference state is known than in the 
case when both reference qubits are in unknown states. 



V. TWO-QUBIT PROGRAM, n COPIES OF 
DATA STATE 

A. Formulation of problem 

We now return to the situation in which we possess 
only one copy of each of the two states we are trying 
to distinguish, but now we have n > 1 copies of the 
unknown state. This means that we want a POVM that 
will distinguish the two n + 2 qubit states 



= IV'^AOlV'i)! 
1*2) = \iPi)a <S> 1^2)1 



\ih) n ®\ih)c, (5.1) 



that occur with a priori probabilities r/i and 772, respec- 
tively. Again we give two protocols, one for minimum 
error discrimination, and one for optimum unambiguous 
discrimination between the states. 

To this end we now define the spaces and operators 
that we will need. Let S be the space of symmetric states 
in H., where Tt is the two-dimensional space for a 

single qubit. £ is an n + 2 dimensional subspace. I^i) 
is an element of E <g) Ti = S\ and is an element of 
Ji (gi £ = 5" 2 . Their intersection, So = S\ n S2, is the 
space of symmetric states in (^) n+2 Tl. Sq is a subspace 
of dimension n + 3. Let K be the subspace of ® n+2 H 
generated by Si and S2. The dimension of K is 3n + 5. 
Let 53 be the orthogonal complement of So in Si, let 
1S4 be the orthogonal complement of So in S 2 , and let 
L be the orthogonal complement of Sq in K . As was 
discussed in Section II, because we do not know what 
and \ip 2 ) are, our problem is to discriminate between the 
density matrices that result from averaging and 
^2) (^2! over l^i) and \ip 2 }- This yields the two density 
matrices 



Pi 



-Py, ® L 



2n + 4 2n + 4 

/'-' = ^-^ p s-2 - 7^—;I®P^ (5.2) 



2n + 4 



2n + 4 



where Ps 1 and Ps 2 are the projections onto Si and S2, 
and and / onto E and Ti., respectively. Consequently, 
we reduced the problem to discriminating between the 
2n + 4 dimensional spaces Si and £2 in K, which is 
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equivalent to discriminating between the n + 1 dimen- 
sional subspaces S3 and £4 in the 2n + 2 dimensional 
space L. 

We will now choose some bases in order to con- 
struct Jordan bases for these subspaces. Jordan bases 
= 0,...A} and {| rj )|j = 0....JV} for two 
N + 1-dimensional subspaces, S p and S r , in general posi- 
tion, are orthonormal bases of their respective subspaces 
({bi)} f° r Sp an( i {l r j)} f° r S r ) that, in addition, satisfy 
( r j\Pk) = $jk cosOk- The angles 9k are called the Jordan 
angles. Now, let |0) and |1) be orthonormal basis vec- 
tors for U. Further, let |uj™ +1) ) (j = 0, . . . n + 1) be the 
unique unit vector in the symmetric subspace of n + 1 
qubits, S, which is the sum of n + 1-tuples with j ones 
and n — j + 1 zeros, 



(n+1) 



(n+1) 



) = I0...0) , 

|0... 01) + |0... 10) + ... + 110...0) 
Vn~TT 



of the above expressions has exactly j ones. Therefore, 
ifj^fc 



(gj\h k ) = 0, 



(5.8) 



and {|5j)} and {|/ij)} form Jordan bases for S3 and S4. 
Let Tj be the two dimensional vector space spanned by 
the nonorthogonal but linearly independent vectors \gj) 
and \hj). The Tj form a decomposition of L into n + 1 
mutually perpendicular two dimensional subspaces. A 
calculation shows that 



, t 1 v n + 1 - j 

{Jj,0\ e 3,0) - n+1 > 

</i-i,i|e;-U> = ^rp 



and 



(ej-i,i|/j.o) = (fj-i,i\ e 3,o) 



(j(n + 2-i)) 1 /2 



n+1 



(5.9) 



(5.10) 



^n+l / 



= 111. 



1) 



(5.3) 



Therefore, 



Then the structure of the two density operators in 
(|5.2|) . in particular the decomposition on the right hand 

side, suggests that we consider \ej. a ) = |itj" +1 ^) ® |a) and 

\fj,a) = I") ® l4" +1) ) where a = 0, 1 and < j < n + 1. 
These vectors form orthonormal bases for Si and S2, re- 
spectively. Let \u^ 1+2 ^) (j = 0, . . . n + 2) be the unique 
unit vector in the symmetric subspace of n + 2 qubits, So, 
which is the sum of n + 2-tuples with exactly j ones and 
11 + 2 — j zeros. This vector can be expressed in terms 
of either the So or Si basis, since it is in both spaces. A 
direct calculation shows that 



( h j\9j) 



n+1 



(5.11) 



and the Jordan angles are all the same. The two den- 
sity operators that we wish to distinguish can now be 
expressed as 



j n+1 

v ' i—l 

j n+1 

P2 = 2(n + 2) [^o + Z>''}vH ( 5 ' 12 ) 
^ ' i—l 



and 



l^" +2) ) ^ J^U\e ifi ) + J-^-Je^), (5.4) 



n + 2 



n + 2' 



\u { r +2) ) = J^^lho) + xl^-M-i,) (5-5) 



n + 2 



n + 2' 



where 



n+2 



^o=Ei4 n+2) )(«r 2)i 
3=0 



is the projection onto Sq. 



(5.13) 



for < j < n + 2. In particular, \u { n+2) ) = |0,0,...,0) 
and |4"+ 2) )) = |1,1,...,1). 
We now introduce the vectors 



and 



\hj) = 



I'//) = \/ ^^l e J'.°) - \l t 2 1 J ' |ej-i,i)» (5-6) 



n + 2 



n + 2 



n + 2 - j 
n + 2 



B. Minimum error discrimination strategy 

For minimum-error identification we have to deter- 
mine the eigenvalues and eigenstates of the operator 
A = r]2P2 — V1P1 which, using Eq. H5.12fl . can be written 
as 



A = ^—[(m-Vl)Ps +V2j2\ h i)( h i 



2n + 4 
-Vi^2\9i){9: 



(5-14) 



for 1 <j < n + 1. The \gj)'s and \hj)'s form orthonormal 
bases for S3 and S4. Each vector on the right-hand sides 



A is diagonal in So and it is straightforward to carry 
out the diagonalization in each of the two-dimensional 
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subspaces spanned by \gi) and \hi), yielding the spectral 
representation, 

n+2 n+1 

A = E A oK)( u jl+E (A+IVt+XVi+l + A-IVi-X^-l) , 

3=0 i=l 

(5.15) 
(5.16) 



where 



A 



7?2 ~ ??1 

2n + 4 



and 



A± = 



1 



4(n + 2) 



('72 - »?i 



/ 2 4n(n+2) 

±4/(772-771) + 



(n + 



^¥«i«2^ .(5.17) 



The eigenvalue associated with Sq, Ao, has a degeneracy 
(n + 3) and the eigenvalues associated with S3 and 5*4, 



X±, have a degeneracy (n + 1) each. Furthermore, in the 
nonorthogonal basis of the two-dimensional subspace Ti 
given by |&) and 



^l + 4-2c±/(n+l)' 



(5.18) 



where 



c± = (n + 1) 



(n+1) 2 ^1 

(2n + 4)A± 



(5.19) 



We find that A_ is unconditionally negative, Ao is nega- 
tive if r] 2 < rji and positive otherwise, and A + is uncon- 
ditionally positive. By making use of Eq. (|2.17l) we find 
that the minimum error probability for identifying the 
state of the data qubits is given by 



J 



Pe — Vmin 



x 2 4n(n+2) 
/ ' (n+1) 2 'Imin'ln 



(5.20) 



where 7] min (Vm ax)js the smaller (larger) of {771, 772}- Ac- 
cording to Eq. H2.16fl . the minimum error probability is 
reached with the help of the optimum detection operators 



n opt = 



J27=l \<Pi-){<Pi 



(5.21) 



Ps a +Y£l\<Pi-)(<Pi-\ if^>| 



which is the projection onto the strictly negative eigen- 
sapce of A, and IT^* = Ik — n° pt , where the identity 

Ik = P Sa + T Ji=i{\ i Pi+)( l Pi+\ + \<Pi+){<Pi+\)- p s is given 
in Eq. (|5.13|) . Clearly, the measurement that identifies 
the state of the data qubits with the smallest possible er- 
ror is a joint projection measurement on all of the qubits. 
It should be noted that for n = 1 the formulas in this Sec- 
tion reduce to those of Sec. III. A whereas for moowe 
have that P E -> rj mm /2 < 1/4. 



C. The optimal universal unambiguous bound 

We now want to consider the unambiguous discrimina- 
tion between the subspaces S% and S 2 in K, or equiva- 
lently between S3 and S4 in L. Let be the orthogonal 
complement of Si in K. is equal to the orthogo- 
nal complement of Si+2 i n L. S t - is an n + 1 dimen- 
sional subspace. The POVM which unambiguously dis- 
tinguishes between Si and S 2 has the form III = aP s ± , 



and n 2 = 0Pg± , where the P's are orthogonal projec- 
tions onto Si or S%, and the a and (3 are positive real 
numbers between zero and one, which are so chosen that 
IIi, n 2 , and n = I — Hi — II 2 are the elements of a 
POVM on K. 

Let us define \g^~) in S^-, and \h^) in S2 by the for- 
mulas 



5i) + — — — \9i ) 



n+1 
1 ,. 



\9i 



n + 1 



n + 1 

on and we have that 



n + 1 



K), (5-22) 



n+1 n+1 

^ - E i^>kH % = E i^x^i- ( 5 - 23 ) 
i=i i=i 

The a and /3 can now be chosen so that Ho restricted to 
each Ti is non-negative. The matrix which represents IIo 
on T i: in the basis {\gi), \gi~)}, is 



1 q » 2 + 2 " V^n+n) 
1 "(n+l) 2 " (n+1) 2 

„ V"("+") 1 a a 

" (n+1) 2 1 P (n+1) 2 



(5.24) 



This matrix must be positive. Therefore, IIi, n 2 , and n 
form a POVM if and only if 
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a-P-ap r +2n 9 >0, 



provided that 



We also have that 



< a,0 < 1. 



Tr(nipi) = 



277 + 2' 



Tr(n 2/ o 2 ) = 



(3n 
2n + 2' 



(5.25) 



(5.26) 



(5.27) 



If we are given pi with probability rji and p 2 with prob- 
ability 772 then the probability that this POVM success- 
fully distinguishes <Si from S2 is 



p — 

± succ — 



2n + 2 



(5.28) 



if the a and /3 satisfy the constraints above. If we set 

l-a 







(5.29) 



(n+l) 2 



in the above formula, which is the maximum allowed by 
(|5.25l) . and differentiate, we find that P SUC c has a maxi- 
mum value when 



r 



a = < 



if m < 



n±L( n + 1 - Jm) if , * 1V , <Vl< TTT^ 
n 2 +n \ y rji J l+(n+l) J — 11 — l+(n+l 



if rji > 



("+ir 

l+(n+l) 2 



The maximum value of P SUC c is 
1 ,n + l 



p — 

* max — 



- {mm) 1 ' 2 ] 



) 2 

(5.30) 
(5.31) 



provided < Vi < l+^+iy^ ■ usin S tne center une 

in (|5.30() . This clearly can be obtained by a POVM only. 
If rji is to the left of this interval (first line in l|5.3Q(| s | 
the optimum measurement is the projection P s ± , which 
unambiguously identifies P2 with a success probability 
Psucc = r]2n/(2n + 2). If 771 is to the right of this inter- 
val (last line in l|5.30l) ) the optimum measurement is the 
projection P s ± which unambiguously identifies pi with 

a success probability Psucc = i]in/(2n + 2). It should 
be noted that the optimum failure probability is given as 
Qp = 1 — Psucc ■ For n = 1 these expressions reproduce 
the corresponding ones in Sec. III.B. For r\\=r\2 = 1/2, 
when their difference is the largest, Qp OVM = 1—P ma x — 



and Qp ^ — 1 Psucc 
00 we have Q p F OVM 



,(1,2) 



1^44, as a function of n. 



2n+4 

For n 00 we have Q F F OVM -> 1/2 and Q ( p' 2) -> 3/4, 
so the POVM outperforms the projective measurements 
quite significantly. Furthermore, Pg < Qp OVM /2 always 
holds, as it should. 

If we use the universal POVM III, II2, and IIo to un- 
ambiguously discriminate between the states \^i) and 



|^2) without averaging over them, we find that P SU cc the 
probability of success is 

Psucc = J7ia<*i|ni|*i)+J72/9(*2|n 2 |* 2 ) 

(ma + mm-\(^M2}\ 2 ), (5.32) 



n + l 



where the a and (3 satisfy the same constraints as above. 
Inserting their optimal values from Eqs. I|5.30|l and (|5.29() 
we find 

~ 1 ~ |^2)| 2 r , , 1/2l q „s 

p o P t = — ; [n + 1 - 2(771772) 7 ], (5.33) 



71 + 2 



with the same restrictions on 771 as in the previous para- 
graph. If we average the overlap term over all possible 
choices of the \ipi)'s we can replace it with its average 
value of 1/2, and we recover (|5.31(l . 

As expected, the optimal success probability is an in- 
creasing function of n. The more copies of the unknown 
qubit we possess, the greater our chance of identifying it. 

The 77, — > 00 limit of P op t can also be achieved by a dif- 
ferent strategy than the one we are employing here. With 
a very large number of copies of the unknown qubit, we 
could employ state reconstruction techniques to find out 
its state |36J. For example, suppose we have determined 
the state of the unknown qubit to be \ipo)- While we know 
what this state is, we do not know if it is equal to the first 
or the second program state, \tpi) or |t/> 2 ) ■ We therefore 
project each of the two program states onto the state or- 
thogonal to the reconstructed state. That is, we take the 
program qubit that we know is in the state |^) (i = 1 or 
2) and measure the projection operator Po± = (ipo\- 
If \ifio) = \ipi), which is given with the a priori probabil- 
ity 77!, then this measurement succeeds (gives 1) with 
a probability of (ipi\ip2)\ 2 , and if \tpo) — Ifo), which is 
given with the a priori probability 772, it succeeds with 
the same probability. Therefore, the total probability of 
success is just (7/1 + r]2)\{"ipi\ip2}\ 2 , which is the same as 
the 77, — + 00 limit of Eq. (|5.33(l . While these strategies give 
the same result for an infinite number of copies, there is 
a difference between them for n finite. The strategy that 
led to Eq. (|5.33|) . will never produce an erroneous result, 
while the strategy based on state reconstruction can. If 
our determination of the state of the unknown qubit is 
not exact, which will, in general, be the case for finite n, 
then \ip^} will not be orthogonal to cither \ipi) or l^), 
and this will lead to errors. 



VI. ONE-QUBIT PROGRAM, n COPIES OF 
DATA STATE 

In this section we return to the case when one of the 
states to be determined is known and the other is un- 
known. However, in contrast to the treatment that was 
presented in Section IV, we are now provided n copies 
of the states to be determined. In other words, we now 
have n data registers, B\, . . . , B ni and in each one we 
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either have a copy of the known state |0) or a copy of 
the unknown state \ip) stored in the program register C . 
Our task is then to decide whether the qubit Bi is in the 
known state |0)s ; , or whether it is in the unknown state 
of the program qubit stored in register C, for all i. 

Thus, we assume that we have a system of n + 1 qubits, 
labeled B\ , . . . , B n , and C where C is the program qubit 
and B\, ... , B n are the data qubits. Qubit C is always 
prepared in the state \tp). Qubits B i} the data qubits, are 
guaranteed to be all prepared in either |0) or all in \ijj), 
but we do not know which of these two alternatives oc- 
curs. The prior probabilities of these two alternatives are 
rji and 772 = 1 — ryi , respectively. Our task is then to find 
the optimal measurement (POVM) that will distinguish 
the two n + 1 qubit states, 



l*i> 
1*2) 



|0)i®. 
1^)1 ® 



|0) n ® |V>)C , 
■\4>)n®\lp) C 



(6.1) 



where we dropped the subscript B for the data registers 
as it leads to no confusion. 

If the state \ip) is completely unknown, we have to find 
the best measurement strategy that is optimal on aver- 
age. Thus, we have to take the average of the input with 
respect to all possible qubit states. The identification 
problem is then equivalent to distinguishing between two 
mixed states, given by the density operators 



p'l 



f>2 



{|*i)(*i|} av 

{|0)®W|^)<0|®W(V»|} , (6.2) 

{|* 2 )<* 2 |} av 

(|^ ( " +1) (^ (n+1) | , (6.3) 

^ J av 



that occur with the prior probabilities 771 and 772, respec- 
tively. The unknown qubit state can be again represented 
using the Bloch parametrization as \ip) — cos(0/2)|O) + 
e 1 ^ sin(0/2)|l), with |0) and |1) denoting an arbitrary set 
of orthonormal basis states. Here 9 and <j> are the polar 
and azimuthal angle on the Bloch sphere. After perform- 
ing the averaging with respect to all possible values of 9 
and <b we arrive at 



p[ = ~|0}®< n >(0|®< B ><8>I. 

J _ ^ psym 

P2 - {n + 2 f Bl ...B nC ■ 



(6.4) 

(6.5) 



where P B B c is the projector to the n + 2 dimensional 
symmetric subspace of the corresponding n + 1 qubits, 
Bi...B n C, and I c = |0)cc(0| + \l)cc(l\- Eqs. 
and (|6.5|l reduce the problem of the programmable state 
discriminator with single-qubit program and n copies of 
the data state to the problem of discriminating between 
these two mixed states. 

At this point it is useful to introduce the following basis 
for the the n + 2-dimensional symmetric subspace of the 



n + 1 qubits 

K) = 



\u n+2 ) = |11 ... 1> . 

We also introduce 

\V2) = 



I0...0) , 

|0... 01) + |0...10) + 



+ |10...0) 



Vn+ 1 



I0...01) 



(6.6) 



(6.7) 



The two nonorthogonal but linearly independent vectors, 
|u2) and life), span a two-dimensional subspace of the 
entire Hilbert space. It will prove useful later on to define 
two other vectors in this subspace as 



and 



\V2) 



U 2 = 



|0... 10) + ... + 110...0) 



1 



VnTT 



v 2 



1 



V2 



(6.8) 



(6.9) 



The two sets, {^2), (^2)} and {|tt2), ^2)}, each form an 
orthonormal basis in the two-dimensional subspace. We 
then have to distinguish between the density operators 
p[ and p' 2 given by Eqs. I|6.4|l and (|6.5|) that refer to the 
first and second alternative, respectively, and that occur 
with the prior probabilities 771 and 772. After reexpressing 
p[ in terms of the basis states \m) (i = 1, . . .n + 2) and 
|«2), defined by Eqs. I|6.6|l and l|6.7|). respectively, the 
density operators to be discriminated read 



P'l = ^OuiXuil + lvaXval): 



P-2 



/n+2 \ 



(6.10) 
(6.11) 



The subsequent treatment proceeds along exactly the 
same lines that we followed in the previous sections. 



A. Minimum-error discrimination 

For identifying the data state with minimum error we 
have to determine the eigenvalues and eigenstates of the 
operator A' = 772 p' 2 ~ ^lPi where p[ and p' 2 are given 
by Eqs. 1)6. 4|) and (|6.5|l . From the explicit expression 
of these two density operators it is clear that A' is di- 
agonal except in the two-dimensional subspace spanned 
by \v 2 ) and |tj 2 ). It is straightforward to carry out the 
diagonalization in this subspace yielding the spectral rep- 
resentation, 



+ 



A' = Xi\ui){ui\ + y Xi\ipi){ipi\ 



i=- 



n+2 

-^K-X^-l, (6.12) 

j=3 



14 



where where 




and 



n + 1 ^ + ± >M - Vi) 2 + ^rV' 2 m 
6.14) c± = — -=■ y - . 

(6.17) 

Here we introduced rj 2 — 2rj 2 /(n + 2) which is the weight 
of p' 2 in the intersection of the supports of the two density 
^ _ V2 /g -^f-N operators to be discriminated. We find that A_ is uncon- 

3 n + 2 ditionally negative, Ai is negative if 772/(11+ 2) < 771/2 

and positive otherwise, and X + ,Xj > for j ' > 3. By 
for J > 3. Furthermore making use of Eq. |2~T7|l we find that the minimum er- 

I \ _ 1-1 ror probability for identifying the state of the data qubits 

W±)=^ V2 \ C±lV2 \ (6.16) is given by 



{'In 



(6.18) 



71+1 



where we introduced r) m i n (r) max ) as the smaller (larger) 
of {771,2772/(71 + 2)}. According to Eq. I|2.16|l . the min- 
imum error probability is reached with the help of the 
detection operators 



III 



opt 



l«i><"i| + \<P-){<P-\ if??i> 



, (6.19) 



n+4 



and IIj 13 * = I — II° pt , where we have to use the identity 



/ = \ui) (ui\ + J2i=3 \ u i) ( U i\+J2t=- Clearly, the 

measurement that identifies the state of the data qubits 
with the smallest possible error is a joint projection mea- 
surement on the qubits Bi (for i = 1, ... ,71) and C. It 
should be noted that for n — 1 the formulas in this Sec- 
tion reduce to those of Sec. IV. A whereas for 71 — * 00 
we have that Pe — > since in this latter case the mixed 
states that we are trying to distinguish become essentially 
orthogonal. The vanishing of the error probability for 
71 — ► 00 is in accordance with the fact that the data state 
can be in principle exactly determined by tomographic 
methods, without any joint measurement, provided that 
an infinite number of copies is available. After the data 
state has been determined, it is of course possible to tell 
without error whether it is equal to the state |0) or not. 



B. Unambiguous discrimination 

Finally we want to determine the minimum failure 
probability for the unambiguous discrimination between 
the states given by (|6.1U|) and (|6.11|l . For this purpose 
we again use the method described in [3(| and 133. Tak- 
ing one of the reduction theorems derived in |26J into 
account, the most general Ansatz for the detection oper- 
ators can be written as 



IIi = a\u 2 )(u2\, n 2 = (3\v 2 )(v2\ + \ui)(ui\, 



n+2 

E 

i=3 



(6.20) 

where 1^2) and \v,2) were given in (|B~5j) and igrHJ), respec- 
tively. Clearly, IIip 2 = Il2p'i = as required for unam- 
biguous discrimination. As follows from Eq. (|2.19|l , these 
detection operators yield the failure probability 



Q'fa 



rjian 



772 n 



772/371 



fail 



2(7i + l) n + 2 (71 + l)(?i + 2) 



(6.21) 



which again has to be minimized under the constraint 
that n = / — IIi — n 2 is a positive operator, in complete 
analogy to our procedure in Sec. IV A. For IIo we obtain 
the expression 



il = \ui)(ui\ - a\u 2 )(u 2 \ + 

-P\V 2 )(V 2 \ + \V2)(V 2 \ ■ 



v 2 )(v 2 \ 



(6.22) 
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The eigenvalues of n are \x \ = 1, \Xj = for 3 < i < n +2, 
and (i± = (2 - P - a ± y/Ja + (3) 2 - 4a(3n/(n + l))/2. 
They all are non- negative provided that f3 < (n + 1)(1 — 
a)/(n+ 1 — na) . In order to minimize Q f a u while keeping 
Ilo a positive operator we therefore choose 



Upon substituting this expression into Eq. i|6.21|l the 
failure probability becomes a function of a alone and it is 
easy to determine its optimum. Taking into account that 
< a < 1, we find that the minimum failure probability 
is obtained when 



OLopt 



^0 



- ■ 



(n+l)(n+2)?7i 



if (71 + 2)771 < 2?7 2 /(n + l) 
if < (n + 2) m < 2(n 4 
if (n+ 2)771 > 2(71 + 1)772 . 



(6.24) 



Substituting these values into H6.23J) yields f3 opt - Using 
a pt and /3 opt in Eqs. (|6.20|) and (|6.22|) yields an explicit 
expression for the optimum detection operators. 

If (n + 2)r)i < 2772/(71 + 1), which implies that 771 < 
2/[2 + (n + l)(n + 2)], we have n° pt = and n° pt = 

1^2) (^2 1 + T^i=z \ u i)( u i\, which means that the optimum 
measurement is a projection measurement on the sub- 
space orthogonal to the span of p[, i. e. a projection on 



its kernel. On the other hand, for (71 + 2)771 > 2(71 + 1)772, 
i.e. 771 > 2(ti+ 1)/(3ti + 4), the optimum measurement is 
a joint projection measurement on the kernels of p[ and 

p' 2 , where II° pt = |u 2 )(u 2 | and n£ pt = Yn=3 \ u i)( u i\- In 
the intermediate parameter region the optimum measure- 
ment is a generalized measurement. The failure proba- 
bility of these optimal measurements can be summarized 
as 



n> — ) 23- + ^2- 

— \ 2 ' n+2 



21 

2 

m "+2. 
'! 2n+2 



(n+l)(n+2) 



^2^2 



if 
if 
if 



< 



2+(n+l)(n+2) 



< 



2+(n+l)(n+2) 
> 2(n+l) 



'/l 



< 



2(n+l) 



2(n+l) + (n+2) 



(6.25) 



'/1 



2(n+l) + (n+2) 



We notice the the above expressions reduce to the cor- 
responding expressions of Sec. IV. A for n — 1, as they 
should. Then, as in that section, it is also true here that 
the benefit of performing the generalized measurement 
is only marginal. To see the closeness of the best PVM 
(projective valued measurement) to the optimal POVM 
we compare their performance in several ways. The two 
PVMs (first and last line in J^HSJ) deliver the same result 
for 771 = 2/(n + 4). In fact, the reduction of the POVM 
failure probability (middle line) compared to those of the 
projective measurements is largest for this value of r/i. In 
Fig. 3 we display the PVM and POVM failure probabil- 
ities for this value of 771 as a function of n. We see that 
the two curves remain close together for all values of n. 
The difference between these two curves as a function of 
n reaches a maximum, however. It is maximal for n = 5 
as displayed in Fig. 4. Finally, in Fig. 5 we display the 
ratio of the PVM failure probability to the POVM fail- 
ure probability as a function of n. Asymptotically, the 
POVM outperforms the PVM by 50%, their ratio tending 
to the limiting value of 1.5. However, as we see from the 



figure, one needs about a 1000 copies of the data state 
to reach the asymptotic region. Since the difference is 
maximal for five copies we can conclude that one does 
not need more than five copies in order to demonstrate 
performance enhancement due to the optimal POVM. 

To close this section we also notice that, in agreement 
with the general relation derived in |35j |. the optimal 
POVM failure probability is always more than twice as 
large as the minimum error probability P' E of the previ- 
ous subsection. 



VII. CONCLUSION 

We have described a number of quantum devices that 
discriminate between two quantum states. We do not 
possess complete information about the states to be dis- 
criminated. Our devices have two inputs, one for the 
qubit whose identity is to be determined, and the other 
for the copies of one or both of the possible states that it 
can be in. In the case that only one of the states is pro- 
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FIG. 3: Comparison of the optimum performances of the 
failure probabilities of the projective measurements (upper 
curve) and the failure probability of the POVM (lower curve) 
vs. the number of copies n for the value of ni = 2/(n + 4), 
when their difference is the largest for the unambiguous dis- 
crimination of one known state from one unknown state when 
n copies of the date state are provided. 
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FIG. 4: The difference between the upper curve and the lower 
curve in Fig. 3 as a function of n. The difference between the 
performance of the PVM and POVM is maximum for n = 5. 



vided, it is assumed that the other state is known, and 
this knowledge is built into the device. The states sent 
into the second input can be regarded as a program. To 
change the set of states between which we are discrimi- 
nating, we do not have to change the device, but merely 
supply it with a different program. 

We want to point out a striking feature of the pro- 
grammable state discriminators in which copies of both 
of the states to be discriminated are provided. Neither 
the optimal detection operators nor the boundaries for 
their region of validity depend on the unknown states. 
Therefore, these devices are universal, they will perform 
optimally for any set of unknown states. Only the proba- 



FIG. 5: The ratio of the upper and lower curve in Fig. 3 as a 
function of n. Asymptotically the failure rate of the PVM is 
50% higher than that of the POVM. 



bility of success for fixed but unknown states will depend 
on the overlap of the states. However, both this expres- 
sion and its average over all possible inputs is optimal. 

The devices described here demonstrate the role played 
by a priori information. All of them have a smaller suc- 
cess probability than one designed for a case in which 
we know both of the input states, and the device for 
two unknown input states has a smaller success proba- 
bility than one designed for the case when we know one 
of the input states. There is a trade off between flexibil- 
ity and success probability. The more of the information 
about the states that is carried by a quantum program, 
the smaller the probability of successfully discriminating 
between the states, but the larger the set of states for 
which the device is useful. This flexibility suggests that 
programmable discriminators will be useful as parts of 
larger devices that produces quantum states that need 
to be identified. 

We conclude our paper by summarizing what we know 
about programmable discriminators with quantum pro- 
grams in which the programs consist of copies of the 
states to be discriminated. The most general problem 
of this type is when we have tia copies of the state of the 
program system A, nc copies of the state of the program 
system C, and tib copies of the state of the data system 
B. In this case, the task is to discriminate two input 
states 



|*l n ) = 



\1>i)a 



(7.1) 



where the subscripts A and C refer to the program regis- 
ters (A contains \xpi) and C contains \ip2}), and the sub- 
script B refers to the data register. Our goal would be 
to optimally distinguish between these inputs, keeping in 
mind that one has no knowledge of and \1fj2} beyond 
their a priori probabilities. The problem in which the 
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numbers of copies of the program states are equal and 
greater than one, but we have only one copy of the data 
state is solved for equal a priori probabilities 19]. The 
problem in which we have only one copy of each program 
state, but an arbitrary number of copies of the data state 
has been solved here. The general problem remains open. 
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